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The possibility of the D-^^i distortion of TiOg octahedra is examined theoreticaUy in 
order to understand the origin of the G-type antiferromagnetism (AFM(G)) and ex- 
perimentally observed puzzling properties of LaTiOs. By utilizing an effective spin and 
pseudospin Hamiltonian with the strong Coulomb repulsion, it is shown that AFM(G) 
state is stabilized through the lift of the t2g-orbital degeneracy accompanied by a tiny 
-Dsd-distortion . The estimated spin-exchange interaction is in agreement with that 
obtained by the neutron scattering. Moreover, the level-splitting energy due to the 
distortion can be considerably larger than the spin-orbit interaction even when the 
distortion becomes smaller than the detectable limit under the available experimental 
resolution. This suggests that the orbital momentum is fully quenched and the rela- 
tivistic spin-orbit interaction is not effective in this system, in agreement with recent 
neutron-scattering experiment. 
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§1. Introduction 

Perovskite-type Ti oxide RTiO^, have recently attracted much interest because of the rich 
magnetic and orbital phases caused by an interplay of spin and orbital degrees of freedom. 
In these compounds, Ti^'*' has a 3(f configuration, and one of the threefold t2g-orbitals is 
occupied at each transition-metal site. The crystal structure is a pseudocubic perovskite 
with GdFeOs-type distortion in which the TiOg octahedra forming the perovskite lattice 
tilt alternatingly. The unit cell contains four TiOe-octahedra, as shown in Fig. 1. With a 
smaller ionic radius of the R ion, the distortion increases with the decrease in the Ti-O-Ti 
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Fig. 1. GdFeOa-type distortion. 



bond angle from 180°. The bond angle can be controlled by use of the solid-solution systems 
Lai_j,Yj,Ti03 or in i?Ti03, by varying the R ions. 

Recently, the magnetic phase diagrams have been investigated intensively in the plane of 
temperature and the magnitude of the distortion. ^"^^ In particular, LaTiOs exhibits a G-type 
AFM (AFM(G)) ground state with magnetic moment of 0.45 jx-ef' and the Neel temperature 
(Tiv) takes about 130 K. With increasing the GdFeOs-type distortion, Tjv decreases rapidly 
and is suppressed to almost zero, subsequently a ferromagnetic (FM) ordering appears. In Y- 
rich systems and in YTiOa in which the GdFeOs-type distortion is relatively large, the system 
shows a FM ground state accompanied by a large Jahn- Teller (JT) distortion. Recently, the 
AFM-FM phase transition with the second-order-like behavior has been understood from 
the strong two-dimensionality in the spin-coupling realized near the transition point. ^'^^ 
However, the origin and nature of AFM(G) state in the end compounds LaTiOs remains 
unresolved. 

Recent model Hartree-Fock study predicts the AFM(G) state with the spin-orbit ground 
state in the absence of static JT distortion.^' However, in contrast to the prediction from the 
spin-orbit ground state, recent neutron-scattering study shows the spin-wave spectrum well 
described by an isotropic spin- 1/2 Heisenberg model with a nearest-neighbor superexchange 
constant J ~ 15.5 meV.^^^ This also suggests the absence of unquenched orbital moments. 
When the orbital degeneracy remains, a FM spin structure with antiferro-orbital ordering 
in which the neighboring orbitals are orthogonal to each other is eventually expected in 
the presence of the transfers and the Coulomb-exchange interaction. Indeed, recent weak 
coupling study shows that in the cubic-crystal field in the absence of the static JT distortion, 
a FM state, out of which two states -^{yz + izx) | and xy t are alternating is favored both 
by the spin-orbit interaction and by the spin-orbital superexchange interaction. However, 
no evidence for the orbital ordering is found in a resonant x-ray scattering study. 

Recently, a possible orbital liquid state was proposed on the basis of small exchange inter- 
action in the orbital sector in the AFM(G) state. However, in a system with strong orbital 
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Fig. 2. D^d distortion. 

fluctuations, the spin and orbital degrees of freedom strongly couple, and both degrees of 
freedom can not be determined independently to each other. In this circumstance, the origin 
of AFM(G) state in LaTiOs is to be clarified in a self-consistent manner. 

In this letter, in order to understand the above properties of AFM(G) state in LaTiOs, 
we propose that the puzzle is solved if we assume the existence of a small Ds^^-crystal field. 
If the D^d distortion exists, TiOg-octahedron is contracted along the threefold direction 
as shown in Fig. 2. As a result, threefold degenerate t2g-leye\s split into non-degenerate 
lower aig-level {(pi) and twofold-degenerate higher e^-levels {ip2, p3). In the distortion with 
[l,l,l]-trigonal axis, the aig representation is -^{xy + yz + zx) and the eg representations are 
■^{2xy — yz — zx) and -^{yz — zx). Since the lowest level has no degeneracy, this distortion 
significantly and sensitively lowers the energy of the electron system in t\g configuration. 
Though the (i-type JT distortion is stabilized in YTiOs or in the compounds with large 
GdFeOs-type distortion due to the effects of covalency between O 2p and unoccupied d- 
orbitals on the i?-site cations, LaTiOa with small GdFeOa-type distortion turns out to be 
unstable to the distortion. The magnitude of the distortion can be denoted by the Dy,^- 
distortion angle u (see Fig. 2). By utilizing the effective spin and pseudospin Hamiltonian 
in the strong Coulomb repulsion, we solve the self-consistent Hartree-Fock equation and 
show that the emergence of the AFM(G) state and the following features in LaTiOs are well 
explained by the effects of D^d distortion even when u is as small as 1°: 

• Tn -130 K, J ~ 15.5 meV, 

• magnetic moment of 0.45 /xb, 

• no detectable JT distortion, 

• no evidence for orbital ordering in the resonant x-ray scattering. 

We analyze why the orbital ordering has not been detected so far, and propose experimental 
ways for testing the validity of our proposal. 
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§2. Formalism 

We start with the multiband d-p model in which the full degeneracies of Ti 3d and O 2p 
orbitals as well as the on-site Coulomb and exchange interactions are taken into account/^ 
The nearest-neighbor d-p and p-p transfers are given in terms of Slater-Koster parameters 
ypd-K, ypda, yppn and Vppa-^^^ The on-site Coulomb interactions are expressed using Kanamori 
parameters, u, u' , j and j' which satisfy the following relations:^^' u — U+^j: u' — u — 2j 
and j — f. Here, U gives the magnitude of the multiplet-averaged d-d Coulomb interaction. 
The charge-transfer energy A, which describes the energy difference between occupied O 2p 
and unoccupied Ti 3d levels, is defined by using U and energies of the bare Ti 3d and O 2p 
orbitals and Sp as A = + U — Sp, since the characteristic unoccupied 3d level energy 
on the singly occupied Ti site is -\- U. The values of the parameters are estimated by 
the cluster-model analyses of valence-band and transition-metal 2p core-level photoemission 
spectra and the analyses of the first-principle band calculations.^'''' We take the values of 
these parameters as A = 5.5 eV, U — 4.0 eV, Vpda — —2.4 eV, V^^^ = 1.3 eV, T^jpo- — 0.52 
eV, VppT^ — —0.11 eV and j — 0.46 eV throughout the present calculation. The effects of 
the GdFeOs-type distortion are considered through the transfer integrals which is defined 
by using the Slater-Koster's parameters^^^ . In the present calculations, we simulate the 
GdFeOa-type distortion by rotating the TiOe octahedra by angle -|-14.1° and —14.1° around 
the [1,1,1] and [-1,-1,1] axes. As a result, the Ti-O-Ti bond angle takes ~ 157°. The effects 
of the distortion are also considered. We integrate out the O 2p orbital degrees of 
freedom in the path-integral formalism to arrive at the effective multiband Hubbard model 
which includes only Ti 3d orbital degrees of freedom. Starting with thus obtained multiband 
Hubbard Hamiltonian, we can derive an effective Hamiltonian in the low-energy region on the 
subspace of states with singly occupied t2g orbitals at each transition-metal site by utilizing 
a second-order perturbation theory. The states of 3d electrons locahzed at the transition- 
metal sites can be represented by two quantum numbers, the z-component of the spin 
and the type of the occupied orbitals. The threefold degeneracy of t2g-orbitals at each site 
can be described by a pseudospin with the spin-1 operator r. We follow an approach similar 
to the Kugel-Khomskii formulation. We express the 3d electron operators in terms of S 
and T to arrive at the effective spin and pseudospin Hamiltonian: i^eff = -f^cry. + ^t2g + H^g- 
The first term iJcry. denotes the effects of Da^-crystal field. This term is obtained from the 
zeroth-order perturbational processes. The second term Hf^^ is obtained from the second- 
order perturbational processes whose intermediate states contain only t2g-orbital degrees of 
freedom. The third term is obtained from the second-order perturbational processes 
whose intermediate states contain e^-orbital degrees of freedom. We neglect the relativistic 
spin-orbit coupling in the model and discuss this issue later. 
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Fig. 3. Energies of each magnetic structures with config. 1 as functions of the I^ad-distortion angle. See 
the text for details. 

§3. Results 

There are four trigonal-axes in the TiOg octahedron, namely, [1,1,1], [1,1,-1], [1,-1^1] and 
[1,-1,-1] with respect to the x, y and z axes. At this stage, there are several possibilities for 
configurations of the trigonal axes at sites 1, 2, 3, and 4. The spin-couphngs between the 
neighboring sites depend on the configuration. There are five configurations which lower the 
AFM(G) spin structure in energy as shown in Table I. In this letter, we examine the AFM(G) 
state with a trigonal-axes configuration in which [l,l,l]-axis is chosen at every sites (config. 
1), and we show that a considerably small Ds^-distortion does strongly stabihze the AFM(G) 
state and the experimentally observed properties are well explained when we assume the 
distortion. It is an interesting problem to study which configuration has the lowest energy, 
which would require an entirely separate work and is left for future discussion. In the D^d 
distortion with config. 1, the lowest level is -^{xy + yz + zx) at every sites. In Fig. 3, we 
plot the energies of several magnetic structures with config. 1 as functions of Ds^-distortion 
angle ui. In the region of a; > 1°, the AFM(G) state is strongly stabilized relative to other 
structures. This may indeed drive the -Ds^ distortion even in the cost of lattice energy. Since 
the lattice elastic-energy is not accounted in Fig. 3, the real energy minimum including the 
elastic-energy cost should be in the range 0° < a; < 8°. In Fig. 4, wc plot the t2g-\cYG[- 
splitting energy due to the Ds^-crystal field (A^jg^) as a function of u). In the region of 
a; > 1°, the value of l^o-^ is sufficiently larger than the coupling constant of the relativistic 
spin-orbit interaction in Ti"^"*" ion (^^ = 0.018 eV),^^-* though this distortion angle uj > 1° 
is substantially smaller than the tilting angle of GdFcOs-type distortion. In fact, it is not 
surprising if such small D^^ distortion has not been detected so far. 

When the level-splitting energy is sufficiently large as compared to the characteristic energy 
of the spin-orbital exchange interaction, the orbital occupation at each site is restricted to 
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Table I. Trigonal-axes configurations which lower the AFM(G) spin structure in energy. 







site 1 


site 2 


site 3 


site 4 




config 
^1 


.1 


[1,1,1] 
■^{xy + yz + zx) 


[1,1,1] 
-^{xy + yz + zx) 


[1,1,1] 
-^{xy + yz + zx) 


[1,1,1] 
-^{xy + yz + zx) 


config 

ip 1 


.2 


[1,-1,1] 
-^{-xy - yz + z:i^ 


[1,-1,1] 
-y=(—xy — yz + zx) 


[1,-1,1] 
-y=(-x'y - yz + zx) 


[1,-1,1] 
-j={—xy — yz + zx) 


config 
^1 


.3 


[1,1,1] 
■^{xy + yz + zx) 


[1,1,1] 

-^{xy + yz + zx) 


[1,1,-1] 
-^{xy-yz-zx) 


[1,1,-1] 

■^{xy-yz-zx) 


config 
^1 


.4 


[1,1,-1] 
■^{xy-yz-zx) 


[1,1,-1] 
-^{xy-yz-zx) 


[1,1,1] 
-^{xy + yz + zx) 


[1,1,1] 
-^{xy + yz + zx) 


config 
^1 


.5 


[1,-1,1] 

■^{-xy - yz + zx) 


[1,-1,1] 
-^{-xy - yz + zx) 


[1,-1,-1] 
-^{-xy + yz + zx) 


[1,-1,-1] 

-^{-xy + yz + zx) 



0.35 
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Fig. 4. Level-splitting energy due to the £>3d-crystal field (A^g^) as a function of the Dsd-distortion angle. 
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Fig. 5. Characteristic transfers which contribute to the spin-exchange interaction. 
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Fig. 6. Spin-exchange interactions estimated under the assumption that orbital occupation is restricted to 
the lowest level ( Jfjcis ^^"^ "^Hois) ^^'^ plotted. Inset shows those calculated using other parameters within 
the error bars of the analyses of photoemission spectra. See the text for details. 



-^{xy + yz + zx) independently of the magnetic structure. At this stage, we can estimate 
the value of the spin-coupling constant in the a6-plane (Jneis) ^^"^ ^^at along the c-axis 
(■^Heis) under this assumption. In Fig. 5, the characteristic transfers which contribute to the 
spin-exchange interaction are presented. Using the notation of transfer integrals as shown 
in Fig. 5, Jg^jg and J^gj^ are expressed in the strong coupling expansion as, 



jab jc A c\ ' ^'6fJ r\ \^ 2 ' ^ 'iJJ ( '\\ 



In Fig. 6, the values of Jgtis ^^^1 Jneis are plotted as functions of u. Both Jg^jg and Jg^jg take 
approximately 0.015-0.022 eV, and they are in agreement with the value obtained by the 
neutron scattering experiment of ~ 15.5 meV. In LaTiOa, the magnetic moment takes 0.45 
/iB, which is reduced from the spin-only moment. Recent optical measurement shows that 
LaTiOs has a considerably small optical gap of ~ 0.1 eV in the vicinity of the metal- insulator 
phase boundary with strong itinerant character. Therefore, in this system, we expect that 
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some amount of orbital and spin fluctuations remain. The reduction of the magnetic moment 
may easily be attributed to these fluctuations. We note that in this context, the calculated 
Jggis and J^gis in the locahzed limit seem rather small at first sight. To elucidate this issue, 
we have examined how the values of J^eis and J^eis change within the error bars of the 
parameters estimated from the analyses of photoemission spectra. In LaTi03 with strong 
itinerant character, we can expect that the Coulomb and exchange interactions are rather 
small. We calculate the values in the case oiU — 3.0 eV, j — 0.50 eV and A = 5.0 eV and 
in the case oi U ^ 3.0 eV, j = 0.40 eV and A = 5.0 eV. The values take ~ 0.03 eV, which 
are larger than the previous estimate (see Inset of Fig. 6). Considering these variations and 
uncertainties, we conclude that the calculated Jgeis and Jggig in the localized limit are in 
qualitative agreement with the experimentally obtained J in the presence of itinerant charge 
fluctuations. 

Here, a question arises: why has the orbital ordering caused by the occupation of the lowest 
level in the D^d distortion not been detected experimentally so far? Detection of the orbital 
ordering by the resonant x-ray scattering is based on the splitting of the Ti Ap levels induced 
by the orbital ordering of the Ti 3d states. When the Ti 3d orbitals are ordered in the 
D^ii distortion, the Coulomb interaction between Ti Ap and occupied 3d orbitals, U^p^sd and 
the interaction between Ti 4p and ligand O 2p orbitals, Uip-2p make opposite contributions. 
Among the threefold Ti Ap orbitals, the orbital directed along the trigonal axis has higher 
energy due to Uipsd. On the contrary, along this axis, C/4p_2p is weakened so that the crystal 
field works to lower the Ti Ap orbital. In addition, the weakend hybridizations between the 
Ti Ap and O 2p orbitals might lower the Ti Ap orbital. This cancellation may lead to the 
absence of detectable orbital ordering in the resonant x-ray scattering in LaTiOs.^^^ The 
orbital ordering accompanied by the distortion would be more easily detected in the 
AFM(G) compounds with relatively large optical gap and small charge fluctuations: PrTiOs 
and NdTiOg. 

In summary, in order to interpret the experimentally observed properties of AFM(G) state 
in LaTiOa, we have examined the effects of possible Da^-crystal field as a candidate for their 
origin. In the presence of the Da^-crystal field, threefold cubic-t2g-levels are split into non- 
degenerate aig- and twofold-degenerate e^-levels. Since the lowest level has no degeneracy, 
this distortion strongly lowers the energy of the electron system in tig configuration. More- 
over, the calculated Dsd-splitting energy A/jg^ suggests that the -Da^-splitting dominates over 
the relativistic spin-orbit interaction in the region of Ds^-distortion angle a; > 1°. The cal- 
culated spin-exchange constant is in agreement with that obtained experimentally. Finally, 
our proposal for the small D^^ distortion should in principle be observed by the displacement 
of O ions if the experimental resolution is sufficient. 
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